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Abstract — We consider a multlcell MIMO uplink channel 
where each base station (BS) is equipped with a large number 
of antennas TV. The BSs are assumed to estimate their channels 
based on pilot sequences sent by the user terminals (UTs). Recent 
work has shown that, as A^ ^^ oo, (1) the simplest form of 
user detection, i.e., the matched filter (MF), becomes optimal, 
(ii) the transmit power per UT can be made arbitrarily small, 
(ill) the system performance is limited by pilot contamination. 
The aim of this paper is to assess to which extent the above 
conclusions hold true for large, but finite N. In particular, we 
derive how many antennas per UT are needed to achieve r; % 
of the ultimate performance. We then study how much can be 
gained through more sophisticated minimum-mean-square-error 
(MMSE) detection and how many more antennas are needed with 
the MF to achieve the same performance. Our analysis relies on 
novel results from random matrix theory which allow us to derive 
tight approximations of achievable rates with a class of linear 
receivers. 

L Introduction 

As wireless networks are inherently limited by their own 
interference, a lot of research focuses on interference reduction 
techniques, such as multiuser MIMO [1], interference align- 
ment |2| or multicell processing [3). Although these techniques 
can lead to considerable performance gains, it is unlikely 
that they will be able to carry the exponentially growing 
wireless data traffic. Due to this reason, a significant network 
densification, i.e., increasing the number of antennas per unit 
area, seems inevitable. One way of densifying the network 
consists in cell-size shrinking, such as the deployment of 
femto or small cells ID, which comes at the cost of additional 
equipment and increased interference. Another, seemingly 
simpler, but also less explored option is the use of very large 
antennas arrays at the base stations (BSs) Q. However, it is 
well known that accurate channel state information (CSI) must 
be acquired to reap the benefits of additional antennas. This 
poses, in particular in fast fading channels, a challenge as the 
number of antennas grows |6|. Thus, massive MIMO is only 
feasible in time-division duplex (TDD) systems where channel 
reciprocity can be exploited. Here, the channels are estimated 
based on orthogonal pilot tones which are sent by the user 
terminals (UTs). In such systems, additional antennas do not 
increase the training overhead and, therefore, "always help" 
Q. Nevertheless, for a given coherence time, the number of 
orthogonal pilot sequences is limited and they must be reused 
in adjacent cells. This leads to channel estimation impairments 
known as "pilot contamination" H]. 



The use of massively many antennas was advocated for 
the first time in ||5l and has since then received growing 
research interest f9l, fTOl, flT|. It is particularly intriguing 
that with an infinite number of antennas, the simplest forms of 
user detection and precoding, i.e., matched filtering (MF) and 
eigenbeamforming, become optimal, the transmit power can be 
made arbitrarily small, and the performance is ultimately lim- 
ited by pilot contamination. Therefore, several works address 
the problem of how the negative impact of pilot contamination 
can be reduced M, M- 

But what is the difference between massive MIMO and 
classical MIMO techniques whose benefits are well-known 
since more than a decade |12| ? This paper tries to give 
an answer to this question. We first provide a definition of 
massive MIMO as an operating condition of cellular systems 
where multiuser interference and noise are small compared to 
pilot contamination. Whether or not this condition is satisfied 
depends in general on several factors: the number of BS 
antennas N, the number of UTs per degree of freedom (DoF) 
offered by the channel (we denote by DoF the rank of the 
antenna correlation matrix which might be smaller than N), 
the signal-to-noise ratio (SNR) and the path loss. We then 
study if more antennas can compensate for suboptimal user 
detection, i.e., how many more antennas does the MF need 
to achieve minimum-mean-square-error (MMSE) detection 
performance. Interestingly, both detectors achieve the same 
performance with an unlimited number of BS antennas. 

The technical contributions of this work are summarized 
as follows: We consider a very general channel model which 
allows us to assign a different path loss and receive antenna 
correlation matrix to each channel from a UT to a BS. As- 
suming a large system limit where the number of BS antennas 
N and the number of UTs per cell K grow infinitely large at 
the same speed, we derive asymptotically tight approximations 
of ergodic achievable rates with a class of linear receivers, 
accounting explicitly for channel training and pilot contamina- 
tion. Simulations suggest that our approximations are valid for 
even small N and K. Our analysis is different from |5 1 which 
assumes that K/N — > 0. However, we obtain their results as a 
special case. Although we do not fully exploit the generality 
of the channel model in the current paper, it enables future 
studies on massive MIMO under more realistic channel models 
with the goal of analyzing the impact of antenna correlation, 
spacing and aperture. Lastly, the techniques developed in this 
work can be directly applied to the downlink llT3l . 



II. System model 



the matrices ^jik are defined as {I < j,l < L,l < k < K): 



Consider a multi-cellular system consisting of i > 1 cells 
with one BS and K UTs in each cell. Each BSs is equipped 
with N antennas, the UTs have a single antenna. The focus of 
this paper is on the uplink without any form of BS cooperation. 
The received baseband signal vector y e C at BS j at a 
given time reads 



B. Achievable rates with linear detection 



(4) 
(5) 



yj 



= Vp^'^jI^i + ^j 



1=1 
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We consider linear single-user detection, where the jth BS 

estimates the symbols Xjm, m = 1, . . . ,K, of the UTs in its 

(1) cell by calculating the inner products between the received 

vector y-j and the linear filters rj.,„ e C . Two particular 

detectors are of practical interest, the matched filter r^|^ and 



where H,, e C is the channel matrix from the UTs ^j^^ ^^^g ^^^^^^^^ ^mmse ^^ich we define as 

in cell / to BS j, x; = [xn ■ ■ ■ xik] ^ CJv (0, Ik) is the 

message-bearing transmit vector from the UTs in cell /, rij ^ 

CAf {0,I]y) is a noise vector and p > denotes the transmit 

SNR. We model the fcth column vector hjik G C of the 

matrix Hji as 



r""" - h ■ 



where A > is a design parameter and 
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(6) 
(7) 



^jik 



R-jifeWjifc 



(2) 
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are deterministic correla- 



where Kjik = rvjifcii,^.;^. 
tion matrices and Wjik ^ CJ\f (0,1^) are fast fading channel 
vectors. The above channel model is very versatile as it allows 
us to assign a different antenna correlation (including path 
loss) to each channel vector Moreover, the matrices Rjik 
do not need to have full rank. This is especially important 
for large antenna arrays with a significant amount of antenna 
correlation due to either insufficient antenna spacing or a lack 
of scattering. 

Remark 2.1: In particular, Q can be used to represent 
a physical channel model with a fixed number of dimen- 
sions or angular bins P as in f^, by letting tijik — 
^ijik [A OjvxAT-p], where A € C and tjik is the 

inverse path loss from UT fc in cell / to BS j. We will use a 



particular form of this model in Section IV 



A. Channel estimation 

During a dedicated training phase (whose length we ignore 
in this work), the UTs in each cell transmit orthogonal pilot 
sequences which allow the BSs to compute estimates Hj^ = 
[hjji • • • , hjj_ff] G C of their local channels Hj%,. The 

same set of orthogonal pilot sequences is reused in every cell 
so that the channel estimate is corrupted by pilot contamination 
from adjacent cells lfT4l . Q. Under these assumptions, BS j 
estimates the channel vector h„,i. based on the observation 



ijk 



y,k 



^jk ■ 



Yl '^j"= 
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-n,k 



(3) 



where Ujk ^ CAf {0,In) and pr > is the effective 
training SNR. Assuming MMSE estimation i fTSl . we can 
decompose the channel as hjjk — hjjfc + ^jjk, where 
hjjk ^ CAf {0,^jjk) is the channel estimate, hjjk ^ 
CAf (OjRjjk — ^jjk) is the independent estimation error and 



E 



HjjH^j + 22 Hj/Hj; 



i^j 



k l^j k 



(8) 



Remark 2.2: Note that a BS could theoretically estimate 
all channel matrices H^; from the observations ^ to further 
improve the performance. Nevertheless, high path loss to 
neighboring cells is likely to render these channel estimates 
unreliable and the potential performance gains are expected to 
be marginal. Our formulation of r^^^ further allows us to 
treat A (and also Zj) as a design parameter which could be 
optimized. A natural choice is A = ^. 

Using a standard bound based on worst-case uncorrected 
additive noise yields the ergodic achievable rate Rjm of UT 
TO in cell j lfT6l : 



i?j„, = E [log2 (1 + 7j>„)] 



(9) 



where the associated "signal-to-interference-plus-noise ratio" 
(SINR) 7j„j is given by ( [TOJ i on the top of the next page. 



III. Asymptotic analysis 

In this section, we present our main technical results. All 
proofs are provided in lfT3l . Under the assumption that N and 
K grow infinitely large while keeping a finite ratio K/N, we 
will derive deterministic approximations 7j,„ of jjm, such that 



7j> 



Ijr 







(11) 



where '— ^' denotes almost sure convergence. We will refer 
to the quantities 7j,„ as deterministic equivalents of 7jm. The 
large system analysis implicitly requires that the channel co- 
herence times scales linearly with K (to allow for a sufficiently 
large number of orthogonal pilot sequences). However, since 
we will use our results as approximations for realistic system 
dimensions, this assumption does not pose any problem. 
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(10) 



We would further like to remark that, in the asymptotic limit, 
the term inside the expectation of the denominator of jjm 
( [TO| l can be arbitrarily closely approximated by a deterministic 
quantity and we do not need to compute the expectation 
explicitly. Moreover, since E[|7jm|] and 7j„i are uniformly 
bounded, the almost sure convergence in ( [TT] i implies by 
IHIJ [Corollary, p. 218]: 



Rj„i - log2 (1 + 7jm) -^ 0. 



(12) 



In the sequel, we assume that the following technical 
conditions hold: 

A 1: limsup^||Rj;fc|| < oo for all j,l,k. 

A 2: liminf^r -^tTRjij. > for all j, I, k. 

Our first result is a deterministic equivalent of the SINR at 
the output of the matched filter: 

Theorem 1 (Matched filter): A deterministic equivalent of 
7jm for the matched filter is given as 
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Our second result is a deterministic equivalent of the SINR 
with MMSE detection: 

Theorem 2 (MMSE detector): A deterministic equivalent 
of "ijm for the MMSE detector is given as 



-MMSE 
JJ7n 



where 



Sin 



tr*«mT;. + i J2Lk t^3lkm + El^j \^, 



pN'^ 



1 



jlm\ 



^,7fe = ^tr*j,fcT, 

fJ-jlkm — T7trR-j7fcT„-, 



2Rc(z9*,,^;,,„)(l + <5,fe)- 1^,7^1'^^-, 



{i + s^kY 



and where 



• Tj = T(A) and dj — [Sji ■ ■ ■ SjkV ^ ^W ^re given by 

Theorem |3] for D = Iat, S = Zj/N, R*. = *jja Vfc, 
. T^ = T'(A) is given by Theorem H for D == I^, S ^ 



T;„ = T'(A), d'^ 



irn- 



— [Sjlwi 

I 



■5: 



jKm 



]"^ = <5'(A) are 



N 



s = Zj/N,e = * 



j]m, 



given by Theorem |4| for D 
Rfc = *jj7c Vfc. 
Theorems [3] and HI can be found in the Appendix. 

Remark 3.1: The theorem can be drastically simplified if 
a less general channel model is considered, e.g. the same 
correlation matrices Hjik for all UTs, no correlation and only 
path loss, Zj = 0. Due to space reasons, we only state the most 
general form here and provide one special case (cf. Corollary 
pb where 7^^^ can be even given closed form. 

Interestingly, the performances of matched filter and MMSE 
detector coincide for an infinite number of antennas: 

Corollary 1 (Performance with infinitely many antennas): 



-MF -MMSE 
'jm ' ' jni 
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where (3jik = lirajy j^tr ^jik, whenever the limit exists. 

Note that 7°^ corresponds also to the asymptotic signal-to- 
interference ratio (SIR) derived in |J5] Eq. (13)]. 

IV. On the massive MIMO effect 

Let us for now ignore the effects of estimation noise, i.e., 
Pt -^ 00, and consider the simple channel model 



H, 



H 



ji 



N 



l^j (13) 



•,NxP 



N columns of 

e (L are 



where A e C is composed of P < 

an arbitrary unitary N x N matrix, Wjj 
standard complex Gaussian matrices and a e (0, 1] is an 
intercell interference factor. Note that this is a special case 
of Q. Under this model, the total energy of the channel 
grows linearly with the number of antennas N and UTs K, 
since E [ti-H^jH^^] ^ ^trAA^^ == KN. The motivation 
behind this channel model is twofold. First, we assume that 
the antenna aperture increases with each additional antenna 
element. Thus, the captured energy increases linearly with N. 
This is in contrast to existing works which assume that more 
and more antenna elements are packed into a fixed volume, 
see e.g. |18|, |19|. An insufficiency of this channel model is 
that the captured energy grows without bounds as iV — > cx). 
However, we believe that linear energy gains can be achieved 
up to very large numbers of antennas if the size of the antenna 
array is scaled accordingly. Second, the number of DoF P 
offered by the channel does not need to be equal to N ifTOl . 
One could either assume P to be large but constant or to scale 
with N, e.g. P = cN, where c £ (0, 1]. In general, P depends 
on the amount of scattering in the channel and, therefore, on 



the radio environment. Under the assumptions made above, 
we obtain the following corollary from Theorem [Tl 

Corollary 2: For the channel model in ( [T3| ) and p^ — > oo, 
7fcm' f^i" ^11 ^' "^' '-^'^ ^^ given in closed form as 

_ \yIB A J- 



L 



(14) 



noise multiuser interference 

where L = 1 + a{L — 1). 

We define the associated rate R^^ as 



a{L-l) 

pilot contamination 



(15) 



One can make several observations from ([14]). Obviously, 
the effective SNR pN increases linearly with N . Thus, if the 
number of antennas is doubled, the transmit power can be 
reduced by a factor of two to achieve the same SNR. Less 
obvious is that the multiuser interference depends mainly on 
the ratio P/K (number of DoF per UT) and not directly 
on the number of antennas. Moreover, noise and multiuser 
interference vanish for iV, P — > oo while pilot contamination 
is the only performance-limiting factor IJS): 

1 



1 



MF 



-^7o 



(16) 



w.p-i-oo, K/JV-fO '" a{L — \) 
We denote by i?oo the ultimately achievable rate, defined as 



^oo = l0g2 (1 + 7oo) = log2 1 



1 



(17) 



a{L-l) 

It is interesting that even with more sophisticated linear 
single-user detection, such as MMSE detection, the ultimate 
performance limit 7^0 cannot be exceeded (see Corollary [T] and 
ifTOl ). Note that without pilot contamination, i.e., for L = 1 or 
a = 0, the SINR grows without bounds as P, iV — > 00. If P 
is fixed but large, the SINR saturates at a smaller value than 
7oo. In this case, adding additional antennas only improves 
the SNR but does not reduce the multiuser interference. Thus, 
with a finite number of DoF, MMSE detection remains also 
for N ^f CO superior to MF. 

Based on the above observations, we believe that it is 
justified to speak about a massive MIMO effect whenever jjm 
is close to 7oo, or in other words, whenever noise and multiuser 
interference are small compared to the pilot contamination. It 
becomes evident from ( [T4| i that the number of antennas needed 
for massive MIMO depends strongly on the system parameters 
P, K, L, a and p. In particular, there is no massive MIMO 
effect without pilot contamination since 700 — > 00. Thus, 
massive MIMO can be seen as a particular operating condition 
in multi-cellular systems where the performance is ultimately 
limited by pilot contamination and the matched filter achieves 
a performance close to this ultimate limit. To make this 
definition more precise, we say that we operate under massive 
MIMO conditions if, for some desired r/ e (0, 1), 



P^P > 77P0 



(18) 



This condition implies that we achieve at least the fraction 77 
of the ultimate performance with the simplest form of single- 
user detection. Replacing R^^ in the last equation by ( [T5] l and 
solving for P/K leads to 



P 
K 



> 



1 a(L-l) 



£M(l + 7oo)''-l] pNL 



Z2 



(19) 



Thus for a given set of parameters p, N, a and L, we can 
find the fraction P/K necessary to satisfy ( fTS] ). For the partic- 
ular channel model considered in this section, also Theorem |2] 
for the MMSE detector can be significantly simplified: 

Corollary 3: For the channel model in ( [T3] l and pr —>■ 00, 
7fcrn^^' foJ" ^11 k,m, can be given in closed- form as 

1 



-MMSE 
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^2-f 



Y ^X 



Z = Aid 



[l + a^iL-l)]{l-2Z) 

Z2 (Z2 - f ) 



P 



1)] 



1-XL- f £2 






XL+^L^) 



^ p 



2[AZ+f (Z2-1)] 

The last result reveals that also the SINR of the MMSE 
detector depends on P, K, N and p only through the ratio P/K 
and the effective SNR pN. Hence, we can use Corollary [3] to 
determine the ratio P/K necessary to satisfy the condition 



R' 



MMSE 



= log (1 + 7' 



MMSEN 



>vRo 



(20) 



which can be very efficiently solved by a simple line search, 
e.g. via bisection. We assume in the sequel A = \/{pN). 

Before we proceed, let us verify the accuracy of the approxi- 
mations R^^ and P^mse f^^ ^^^^^ (^^ j.^-, jj^ pjg jJJ ^g ^jgpj^^ 

the ergodic achievable rate Rjm of an arbitrary UT with MF 
and MMSE detection as a function of the number of antennas 
N for K = 10 UTs, L = 4 cells, p = OdB and intercell 
interference factor a = 0.1. We compare two different cases: 
P = N and P = N/3. As expected, the performance in the 
latter scenario is inferior due to stronger multiuser interference. 
Most importantly, our closed-form approximations are almost 
indistinguishable from the simulation results over the entire 
range of N. 

Figs. [2] and [3] show the necessary DoF per UT P/K for 
a given effective SNR pN to achieve a spectral efficiency 
of rjRoo with either the MF (solid lines) or MMSE detector 
(dashed lines). We consider L = 4 cells and an intercell 
interference factor a — 0.3 and a ~ 0.1, respectively. The 
plots must be understood in the following way: Each curve 
corresponds to a particular value of ry. In the region above 
each curve, the condition ([T8]l, respectively (pOJi, is satisfied. 
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odic achievable rate with MF and MMSE detection versus number 
N for P = {N, N/3} for p = dB. 
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Fig. 2. Degrees of freedom per UT P/K necessary to achieve »yiJoo versus 
effective SNR pN for L = 4 and a = 0.3. 



Let us first focus on Fig. l2]with a — 0.3. For an effective 
SNR pN = 20 dB (e.g. p = OdB and TV == 100 == 20 dB), 
we need about P/K — 90 DoF per UT to achieve 90% of 
the ultimate performance Roo, i.e., 0.9 x 2.2 w 2b/s/Hz. If 
P Ki N, only a single UT could be served (Note that this 
is a simplifying example. Our analysis assumes K :^ 1.). 
However, if we had N = 1000 = 30 dB antennas, the 
transmit power p could be decreased by 10 dB and 10 UTs 
could be served with the same performance. At the same 
operating point, the MMSE detector requires only ~ 60 DoF 
per UT to achieve 90% of the ultimate performance. Thus, 
MMSE detection would allow us to increase the number of 
simultaneously served UTs by a factor ||| = 1.5. This example 
also demonstrates the importance of the relation between N 
and P. In particular, if P saturates for some N, adding more 
antennas increases the effective SNR but does not reduce the 
multiuser interference. Thus, the number of UTs which can be 
simultaneously supported depends significantly on the radio 
environment. We can further see that adding antennas shows 
diminishing returns. This is because the distances between 
the curves for different values of 77 grow exponentially fast. 
Remember that for 77 = 1, a ratio of P/K = cx) would be 
needed. A last observation we can make is that the absolute 
difference between MF and MMSE detection is marginal for 
small values of t] but gets quickly pronounced as ry ^- 1. 

Turning to Fig. [5] for a — 0.1, we can see that for the 
same effective SNR pN — 20 dB and the same number 
of DoF per UT P/K ~ 90 as in the previous example, 
only 80 % of the ultimate performance are achieved by the 
MF. However, since the intercell interference is significantly 
smaller compared to the previous example, this corresponds to 
0.9x5.1 w 4.6b/s/Hz. Thus, although we operate further away 
from the ultimate performance limit, the resulting spectral 
efficiency is still higher. With MMSE detection, only 35 
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Fig. 3. Degrees of freedom per UT P/K necessary to achieve TyiJoo versus 
effective SNR pAf for L = 4 and « = 0.1. 



DoF per UT are necessary to achieve the same performance 
and, consequently, 90/35 « 2.5 times more UTs could be 
simultaneously served. With decreasing intercell interference 
(and hence decreasing pilot contamination) the advantages of 
MMSE detection become more and more important. 

V. Conclusions 

We have studied the uplink of a MIMO multi-cell system 
where channel estimates are prone to pilot contamination. 
Under a very general channel model, allowing for an individual 
correlation matrix for each UT, and the assumption that the 
number of BS antennas and UTs grows large, we have derived 
deterministic equivalents of achievable rates with matched 



filter and MMSE detector. These results have been used to 
study the performance of both detectors in the large but finite 
N regime. Interestingly, their performance depends mainly on 
the DoF per UT the channel offers and the effective SNR. 
Moreover, we have determined (i) how many antennas are 
needed to achieve 77 % of the ultimate performance limit and 
(ii) how many more antennas are needed by the MF to achieve 
MMSE performance. Our results can be also applied for the 
study of the downlink 1 13| and could be further used to analyze 
the effects of more realistic channels models, such as antenna 
correlation, spacing and aperture. 

Appendix 

Theorem 3 (JM Theorem ]]): Let D G C^""^ and S e 
C be Hermitian nonnegative definite matrices and let 

H € C be a random matrix with columns h^ = 

-j^R^Uk, where u^ G C are random vectors of i.i.d. 
elements with zero mean, unit variance and finite 8th order 
moment, and R^ S C are deterministic covariance 

matrices. Assume that D, S and R^,, k = 1,...,K, have 
uniformly bounded spectral norms (with respect to N). Let 
N,K ^ 00, such that < liminf ^ < limsup^ < 00. 



Then, for any p > 0, 

^trD(HH^ + S + plAr)"' 
where T(p) e C is defined as 
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trDT(p)-^0 



T(P) 
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Rfc 



N ^-^ 1 
fe=i 



hip) 



S + pI 



N 



and the following set of K implicit equations 
1 



Skip) 



N 



trRfcT(p), fc = l,. 



,K 



has a unique solution dip) — [Siip) ■ ■ ■ Sxip)] > 0. 

Theorem 4: Let G C be a Hermitian nonnegative 

definite matrix with uniformly bounded spectral norm (with 
respect to N). Under the same conditions as in Theorem l3] 
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N 



ti-D(HH" + S + pIjv) ^©(HH^^ + S + pIjv)" 



N 



ti-DT'(p)^0 



where T'ip) G C is defined as 



T'ip)=Tip)&Tip) 



<ta 



^kS'.ip) 



rT(p) 



with T(p) and S^ip) as defined in Theorem p] and S'ip) 
[S'iip)---SKip)V given by 

S'ip)^ilK-Jip))-'vip) 
itrRfcT(p)RiT(p) 



[Jip)]ki = 



Nil + Skip))' 



[v(p)], = -trRfeT(p)eT(p) 
where J(p) G C and v(p) G C . 
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